Dynamics of a stochastic SEIQR model driven by Lévy jumps with bilinear incidence rates

In this study, we propose a stochastic SEIQR infectious disease model driven by Lévy noise. Firstly, we study the existence and uniqueness of the global positive solution of the model by using the stop-time. Secondly, the asymptotic behavior of the stochastic system at disease-free equilibrium and endemic equilibrium are discussed. Then, the sufficient condition for persistence under the time mean is studied. Finally, our theoretical results are verified by numerical simulation.


Introduction
Infectious diseases have always been one of the important threats to human health, and the control of infectious diseases is an important issue in human society.It is well known that Kermack and McKendrick first proposed the SIR Model based on the Indian plague model [1].Piovella [2] proposed a SEIR model considering the type E(t) that receives but does not propagate.
Most of the previous models of infectious diseases were basically considered on the basis of assuming the free movement of individuals in the population, and rarely considered the problem of having isolation chambers.With the onset of COVID-19 in 2020, the prevention and treatment of infectious diseases has become one of the topics of research for governments around the world.After the continuous exploration of prevention, the introduction of isolation chamber Q has an excellent effect on timely controlling of infectious diseases.Therefore compared with the previous SIR, SIRS, SEIR and other models, SEIQR model can more accurately describe the prevention and control of infectious diseases.Liu et al. [3] proposed a deterministic SEIQR(S: Susceptible; E: Exposed; I: Infected; Q: Quarantined; R: Removed) epidemic model: where t is the time; and the normal number A represents the number of births and immigrants entering S(t) per unit time; α is the proportion from S(t) to E(t); b is the incidence of incubation period; μ, a, d, k represents removal rates from classes S, E, I, Q respectively; h stands for isolation strength; δ is category I natural mortality; n is the natural mortality rate of category R recoveries due to immune system impairment, age and other reasons.
Converting model (1) to the following form: , where 1 � i, j � m, FV −1 is called a regenerative matrix.The basic regeneration number is the spectral radius of the regeneration matrix.Basic reproduction number R 0 of system (1) is Since the infectious disease model is affected by many unpredictable environmental noises, adding random interference to the deterministic model can reflect the transmission law more accurately.In [4][5][6][7][8][9][10], white noise interference factor was added to the deterministic model to study the dynamic behavior of a stochastic infectious disease model.Gaussian white noise is used to describe a class of relatively stable, continuous random interference.Tornatore et al. [4] proposed a stochastic SIR model with or without distributed time delay and studied the stability of disease-free equilibrium.Xu et al. [5] studied a kind of SIRS model, proved the existence and uniqueness of the positive solution of the model and obtained the conditions of disease extinction for epidemics.Zhao [6] studied the relationship between the threshold value of stochastic SIRS model with saturation incidence and the extinction and persistence of epidemic diseases.Hieu [7] mainly studied the stochastic SIRS model under telegraph noise and gave the conditions of disease persistence and disease-free equilibrium stability.Cai [8] mainly discussed the limit of transforming SDE model to discrete-time system and proved that the regeneration number can be used to judge the relevant properties of SDE model by using Markov semigroup theory.Yuguo et al. [9] analyzed that the distribution of stochastic SIR model solutions is absolutely continuous.Liu et al. [10] demonstrated that the system has a unique global positive solution and established sufficient conditions for disease persistence.Hattaf et al. [11] proposed and analyzed a stochastic SIR Epidemic model with specific functional response and time delay, and compared the difference of the basic regeneration number between the deterministic model and the stochastic model.Similarly, Lan et al. [12] studied a stochastic SIS model with saturated exposure rates and also found that the conditions for extinction of the disease were much weaker than the corresponding deterministic model.Ali and Khan [13,14] studied the dynamic properties of stochastic SEIR and SIRS models with saturation rate and simulated them using Legendre spectrum method.
However, disease can be affected by a variety of natural mutations, such as volcanic eruptions, chemical pollutants, and sudden climate changes, which are often not accurately described by stochastic models of Brownian motion.Therefore, many studies on natural mutation factors will use Le ´vy jump to describe.This perturbation can more accurately describe the impact of mutation factors, and more deeply understand and predict the trend of disease spread and development.According to the Le ´vy-Ito ˆdecomposition theorem [15], Le ´vy noise is composed of Brownian motion, independent Poisson random measures, and deterministic drift terms, so Le ´vy noise has a wider applicability than white noise [16][17][18][19][20]. Zhang and Wang [16,17] studied SEIR model and S-DI-R model driven by white noise and Le ´vy noise respectively.Zhang et al. [18] studied the dynamics of a stochastic SIS epidemic model with saturation incidence and dual epidemics and obtained sufficient conditions for the average extinction and persistence of both epidemics.Liu et al. [19] discussed the persistence and extinction of a delayed vaccination SIR epidemic model with temporary immunity and Le ´vy jumps and analyzed the relationship with the intensity of Le ´vy noise and the duration of vaccination with the duration of disease and the duration of eradication.EL Koufi et al. [20] considered a stochastic SIR model with a saturated incidence rate and saturated treatment function incorporating Le ´vy noise.Based on this, a stochastic SEIQR model with Le ´vy jumps can be considered: where B i (t) represents standard Brownian motion with filter {F t } t>0 on a complete probability space (O, F, P), and they are independent of each other; where

Existence and uniqueness of the global positive solution
We assume that the jump diffusion coefficient satisfies the following conditions: (H1): R Lemma 2 Suppose that Conditions (H1) and (H2) hold, for any given initial value ðSð0Þ; Eð0Þ; Ið0Þ; Qð0Þ; Rð0ÞÞ 2 R 5 þ , system (2) has a unique positive solution (S(t), EðtÞ; IðtÞ; QðtÞ; RðtÞÞ 2 R 5 þ , and the solution lies in R 5 þ with probability 1. proof According to (H1), system (2) satisfies the local Lipschitz condition, for any given initial value, the system (2) has a unique local solution (S(t), where τ e represents the blasting time.To prove the existence of a global solution, just prove τ e = 1, a.s.
Let k 0 be a sufficiently large positive number such that the initial value (S(0), E(0), I(0), Q , For any k � k 0 , the stopping time is defined as: Obviously, τ k is monotonically increasing with respect to k and If we could prove that τ 1 = 1, then τ e = 1.
Next, we use the proof by contradiction to prove.Suppose that hence there exist constants N > 0 and ε 2 (0, 1) such that P(τ 1 � N) > ε.
Then there exists an integer where n is a normal number to be determined later.It is easy to judge that u + 1 − ln u > 0 is true, then V > 0. And according to Ito ˆformula, we get: where System ( 2) is substituted into the above formula and the fundamental inequality is used, then According to Taylor's formula and (H2),  Thus τ 1 = 1 a.s., then (S(t), E(t), I(t), Q(t), R(t)) will not explode in a finite amount of time.Therefore, there is a globally unique positive solution for system (2).
Asymptotic behavior around disease-free equilibrium of the deterministic model is the disease-free equilibrium point of the deterministic model.When R 0 � 1, P 0 is asymptotically stable.Next we will discuss the asymptotic behavior of the solution of the stochastic model at the disease-free equilibrium point.Theorem 3 Suppose that conditions (H1), (H2) hold, if R 0 � 1, and the following conditions are met: For any given initial value ðSð0Þ; Eð0Þ; Ið0Þ; Qð0Þ; Rð0ÞÞ 2 R 5 þ , the solution of system (2) has the following properties: where proof Define the following functions: Thus Define the Lyapunov function again: According to the Ito ˆformula, then where For the Lyapunov function to be asymptotically stable, then m 1 , m 2 , m 3 , m 4 , m 5 > 0. Thus m 1 ; m 2 ; m 3 ; m 4 ; m 5 ; m; M are described in Theorem 3. Integrating both sides of (10) from 0 to t and taking the expectation: This completes the proof of Theorem 3.

Remark 1:
Theorem 3 shows that, with some suitable conditions, the solution of system (2) fluctuates around the disease-free equilibrium point P 0 , the wave intensity is related to noise intensity σ i and C i .The bigger the σ i and C i , the bigger the fluctuation.That is, the greater the random disturbance, the farther away the solution of system ( 2) is from the disease-free equilibrium point P 0 of the deterministic model, at which time the disease will disappear.Next, we will verify the correctness of Theorem 3 through numerical analysis(see From the observation of Fig 1, it can be seen that under certain conditions of parameters, the system will stabilize in a situation where only susceptible persons exist.The infected person, exposed person of the virus will disappear, which means the disease will disappear.

Asymptotic behavior around endemic equilibrium of the deterministic model
P* is the endemic equilibrium point of the deterministic model.When R 0 > 1, P* is asymptotically stable.Next we will discuss the asymptotic behavior of the solution of the stochastic model at the endemic equilibrium point.Theorem 4 Suppose that conditions (H1), (H2) hold, if R 0 > 1, and the following conditions are satisfied: For any given initial value ðSð0Þ; Eð0Þ; Ið0Þ; Qð0Þ; Rð0ÞÞ 2 R 5 þ , The solution of system (2) has the following properties: where proof Define the following functions: According to the Ito ˆformula, then where For the Lyapunov function to be asymptotically stable, then l 1 , l 2 , l 3 , l 4 , l 5 > 0.Where l 1 , l 2 , l 3 , l 4 , l 5 , L are described in Theorem 4. Integrating both sides of (13) where l is described in Theorem 4. This theorem is proved.Remark 2: Theorem 4 shows that when some conditions hold, the solution of system (2) oscillates around P*, and the intensity of the vibration is related to the noise intensity.When the degree of disturbance is greater, the solution of the system (2) is further away from the local equilibrium point of the deterministic model, and the disease will persist.
Next, we will verify the correctness of Theorem 4 through numerical analysis(see

Persistence of disease
Theorem 7 Assume that S(t), E(t), I(t), Q(t), R(t) is the system (2) with initial value (S(0), E (0), Integrating system (2) from 0 to t, dividing by t and substituting F(t) yields: According to the model (2),  20) into ( 16): According to FðtÞ ¼ < SðtÞ >¼ 1 t where ε is any positive constant.In combination with (24) and (25), This theorem is proved.Theorem 7 states that under certain conditions, the disease will continue to spread.This means that the disease persists among the population and is not conducive to further management.